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Abstract 

The non-perturbative input necessary for the determination of the 0{g^) part of 
the weak coupling expansion of the free energy density for SU(2) and SU(3) gauge 
theories is estimated. Although the perturbative information completing the con- 
tribution to this order is missing, we give arguments that the magnetic fluctuations 
are dominated by screened elementary magnetic gluons. 



1 Introduction 



Perturbation theory in its original formulation fails beyond 0{g^) in the calculation 
of the free energy density of non-Abelian gauge theories 0. Very recently Braaten 
has pointed out that the coefficients of the higher powers in the weak coupling 
expansion can be determined by invoking non-perturbative information from the 
effective theory of three-dimensional static magnetic fluctuations. 

In Ref. a systematic two-step separation of the perturbatively treatable fluc- 
tuations from the static magnetic sector has been proposed. In the first step the full 
(electric and magnetic) static sector is represented by an effective three-dimensional 
theory. In this theory a massive (m^;) adjoint scalar field stands for the screened 
electric fiuctuations. In the second step this theory is matched onto an effective 
magnetic theory (MQCD) with a separation cut-off Am- While the contribution of 
the non-static and of the static electric modes to the free energy can be safely calcu- 
lated perturbatively (expansion parameters are g(T) and itie/T, respectively), the 
magnetic sector is still non-perturbative and should be investigated by numerical 
methods. 

The high temperature free energy of the full theory can then be written in addi- 
tive form: 

/qcd(T » T,) = fNsA^M) + /M(AAf), gT > Am > g'T. (1.1) 

The double subscript NS, E refers to the fact that the first term contains contri- 
butions to the free energy from the non-static as well as from static electric type 
fiuctuations. 

The strategy of the hierarchical calculation has been tested by Braaten and Nieto 
by reproducing from the effective electric QCD theory the 0{g^) term in the free 
energy density of the SU(N) gauge theory, which earlier has been calculated directly 
in the full theory by Zhai and Kastening ||^. This is fully part of the term fNS,E- 

The necessary non-perturbative information to go beyond this stage comes from 
the minimal MQCD theory. The purpose of our paper is to present first results 
of a lattice analysis of MQCD, ie. the 3-dimensional SU(N) gauge theory. This 
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does provide the non-perturbative input needed to evaluate the 0{g^) part of /m, 
denoted in the following by 



(0) 
M ■ 



This piece of information can be derived from the simplest three-dimensional 
gauge theory: 
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where 



r(0) 

^MQCD 



[1.3) 



is the only operator of dimension four contributing to leading order. The gauge 
coupling of MQCD with the necessary accuracy is given by = g^T. It is worth to 
note that the leading {0{g^)) correction to g^ comes from the finite wave function 
renormalization factor due to the integration over the static electric field. 

The determination of f[f represents a certain interest in itself. It provides 
information about the nature of quasi-particle excitations in the high temperature 
phase. Similar to the situation in the electric sector of the SU(N) gauge theories 
one also may expect that at very high temperature a weakly interacting gas of 
some (quasi)particles dominates the free energy density contribution of the magnetic 
sector. Yet, the nature of the "constituents" of this gas is still to be clarified. 

Inspired by the reduction strategy we have recently investigated the problem 
of magnetic screening both in pure SU(2) and in SU(2)-Higgs models In Lan- 
dau gauge we have analyzed the vector (Ai) two point correlations and determined 
the propagator mass from the corresponding Euclidean propagator. The numerical 
results were found to be similar to those obtained in analytical calculations from 
coupled gap-equations for the Higgs and the vector channels 0. In the pure gauge 
theory and in the SU{2) gauge-Higgs systems also heavier excitations have been 
identified numerically [§, ||, |lOl and interpreted analytically ||lT| . A unified interpre- 
tation of these rich spectra is, however, missing at present. In particular one has 
to clarify whether the light excitations found so far only in gauge dependent corre- 
lation functions really will dominate the thermodynamics in the high temperature 
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ideal gas limit. Understanding the nature of (quasi)particle constituents in the high 
temperature limit still is one of the major challenges in non-abelian gauge theories. 
The single number offers an interesting input to this discussion as it is sensitive 
to the mass of the thermodynamically relevant lightest magnetic excitations. 



2 The 3-d magnetic free energy density 



Simple analysis of the vacuum diagrams to 4 loop order in the 3-d SU(N) gauge 
theory with cut-off Am leads to the following dependence of /|°^ on the magnetic 
separation scale: 
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(2.1) 



It is worth to remark, that the appearence of the gauge coupling g\ under the 
logarithms implicitly assumes the existence of an infrared scale proportional to it. 

The complete free energy density does not depend on Aj\/. Therefore the coef- 
ficients of the terms diverging with Km should coincide with their perturbatively 
calculable values. Only in this way the Am dependence can cancel against the fully 
perturbatively determined fNS,E- The finite part proportional to 04, however, is 
fully non-perturbative, since all loop diagrams at and beyond four loops contribute 
to it. This is the quantity we would like to extract in a non-perturbative calculation 
on the lattice. 

Starting from 0{g^), contributions to the free energy density are influenced also 
by the higher dimensional operators appearing in the effective magnetic theory. This 
implies, that the 0{g'^) nonperturbative contribution reflects the 0{g^) correction 
of the g1 = g'^T relation and is still calculable within the minimal effective theory. 

The actual procedure for the determination of the coefficients in (|2.1| ) amounts 
to measuring the coefficients of the weak coupling expansion of the internal energy 
density of the system. Since the temperature dependence appears in this theory 
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exclusively through g^, one finds for the energy density 



V dT dgl 

- -r^^.. (2.2) 



The structure of the 3- dimensional energy density and its cut-off dependence is given 
by 

es = + a2h\i + 2a^kMgl + (3a4 - a'^)gt + ^a'^gt In ^ . (2.3) 



The lattice regularization of the minimal 3-d SU(N) gauge theory is defined in 
the standard way: 

Slmqcd = /53 E (l - ^Tr(f/p + f/t)) , /53 = (2.4) 



where Up denotes the Wilson plaquette variable defined in terms of SU(N) valued 



variables Ux,i IT2[- The partition function is given by 



Zlmqcd = / n ^^^^i exp{- Slmqcd)- (2.5) 



The internal energy and the plaquette expectation value 



W = (l-5^T.(C/p + C/^)> = -iA.nZ (2.6) 



can be simply related: 



63 = -3Ai,^(P). (2.7) 

93 



Here A^/ = a ^ is chosen to coincide with the cut-off of the lattice regularized theory. 



3 Plaquette Expectation Value 



The basic non-perturbative input from a lattice calculation is obtained through 
an evaluation of the plaquette expectation value for a SU(N) gauge theory. Its 
perturbative expansion for large j3 has been calculated up to 0{(3~'^) for arbitrary 
dimensions d and on finite lattices [13].0 



ECn,d/3-" (3.1) 
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where Id = — and is a numerical coefficient, which has a weak volume 
dependence. On an infinite lattice a direct evaluation gives = —0.000103 for d = 4 
and as = —0.00095 for d = 3. We note that at order the dominant contribution 
to the expansion coefficients comes from diagrams which are proportional to c?""". 
This allows to estimate also the expansion coefficient at 0{j3^^) which so far has 



only been evaluated in four dimensions |14|. These coefficients are 03,4 = 0.143055 



in the case of SU(2) and 2.960467 in the case of SU(3). Multiplying with a factor 
(4/3)^ one finds as an estimate in three dimensions 

rO.34, SU(2) 
"^''=^-17.02, SU(3) • 

The possible appearance of a logarithmic /9-dependence in the next order expansion 
coefficient, 04^3, follows from combining ( |2.3D and (|2.7|). Its consequences will be 
discussed below. 

As we are finally only interested in the finite part in our aim is to extract the 
coefficient 04,3 in the expansion of (P) . In order to do so we subtract the known part 
of the perturbative expansion from the numerical results for the plaquette expecta- 
tion values obtained from a Monte Carlo simulation and determine the coefficients 
C4 3 and C3 3 from a fit to these differences. We have calculated plaquette expectation 



^We use f3 for the coupling in arbitrary dimensions. 
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SU(2) 



/33 


(P) 


A 


# iterations 


6.0 


0.1752315(78) 


0.4593(17) 


70000 


6.5 


0.1609606(79) 


0.4475(22) 


60000 


7.0 


0.1488726(47) 


0.4413(16) 


50000 


7.5 


0.1384806(60) 


0.4338(25) 


80000 


8.0 


0.1294547(50) 


0.4276(26) 


100000 


9.0 


0.1145530(39) 


0.4248(29) 


110000 


10.0 


0.1027277(40) 


0.4121(40) 


90000 


11.0 


0.0931350(42) 


0.4164(56) 


110000 


12.0 


0.0851753(35) 


0.4060(61) 


100000 


13.0 


0.0784784(42) 


0.4094(93) 


50000 


14.0 


0.0727477(30) 


0.3817(83) 


50000 



Table 1: Plaquette expectation values and normalized differences as defined in 
Eq. (|H73|) calculated on a 16^ x 64 lattice. 



values at a large set of /^a values both for the SU(2) and SU(3) gauge theory, using 
lattices of sizes 16^ x 64 and 32^, respectively . Results for (P) and the difference 



(3.3) 



are summarized in Tables |I| and ^ The difference A is shown in Figure ^ In A 
the leading powerlike divergences have been eliminated. These would anyhow be 
canceled in the final physical result for the energy density through a corresponding 
perturbative calculation within the effective electric theory (NS,E). 

The remaining cut-off dependence results from a possible logarithmic cut-off 
dependence of the free energy density, ie. the term proportional to 04 in Eq. (p.3|). 
We thus expect A to depend on /?3 as follows: 



A = C3,3 + U In Ps + c,,s) + ^ + 0(A 



(3.4) 



The coefficient C4 is directly related to a'^ and could in principle be fixed through a 
perturbative calculation within the effective electric theory. We also note that 03^3 



can directly be evaluated in lattice perturbation theory (see Eq. p7^) 
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1 nnnn 
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15 


0.1892309(59) 


9.423(20) 


10000 


16 


0.1765317(56) 


9.228(23) 


10000 


17 


0.1654524(52) 


9.074(26) 


10000 


18 


0.1556935(48) 


8.931(28) 


10000 


19 


0.1470324(32) 


8.808(22) 


20000 


20 


0.1392884(31) 


8.673(25) 


20000 


21 


0.1323326(28) 


8.618(26) 


20100 


22 


0.1260367(27) 


8.511(29) 


20000 


23 


0.1203172(21) 


8.426(26) 


30300 


24 


0.1150972(20) 


8.351(28) 


30500 


26 


0.1059145(23) 


8.242(40) 


20000 


28 


0.0980936(21) 


8.148(46) 


19000 


30 


0.0913523(15) 


8.085(40) 


36000 



Table 2: Plaquette expectation values and normalized differences as defined in 
Eq. ( p.3| ) calculated on a 32^ lattice. 
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Figure 1: The subtracted plaquette expectation values A as defined in Eq. (|3.3|) for 
tlie SU(2) (a) and SU(3) (b) gauge tlieory. Tlie solid curves show fits with c'^ = 
(case I) and the dashed line corresponds to C4 7^ (case II). For the SU(3) case both 
fits coincide in the coupling range shown. 

Our numerical data for the plaquette expectation values are limited to a /5-range 
that varies by a factor of about 2.5. This is not large enough to be sensitive to the 
logarithmic term given above and at the same time allow a control over the sub- 
leading corrections proportional to C5 3. We thus have analyzed the numerical data 
using two different fits. Since for the moment the 4-loop perturbative information 
on the coefficient of the logarithmic term as well as a direct evaluation of C3 3 are 
missing we had to fit these. In order to see the sensitivity of 04,3 to the variation of 
C4 we have investigated also another option, corresponding to an ideal decoupling of 
the magnetic fluctuations from the rest of the system: C4 = 0. It will become clear 
in the next section that the range of C4 values obtained in this way puts already 
a rather stringent bound on the mass of magnetic excitations in the plasma phase. 
We have fitted A with the ansatze 

_ f C3,3 + 1^ + If case I , 

^ ~ I C3,3 + (c'^ In p3 + 04,3) + ^ case II . ^^'^^ 

Results of these fits to the plaquette expectation values are summarized in Table |^. 
As our data for the SU(2) gauge theory did not show any significant curvature in 
the couphng range explored by us we have fixed 05,3 to be zero in both fits. 
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SU(2) 


SU(3) 




case I 


case II 


case I 


case II 


C3,3 


0.351(5) 


0.45(5) 


7.30(11) 


8.11(20) 


C4,3 


0.635(37) 


1.4(4) 


15.2(3.6) 


99(6) 


c'. 




-0.75(35) 




-29.2(3.4) 


C5,3 






252(29) 





Table 3: Results of the fits using the two fitting functions defined in Eq. (p.5|) . 

We note that the expansion coefficients for SU(2) gauge group are substantially 
smaller than in the case of SU(3). This is in agreement with the expected A^- 
dependence for the expansion coefficients for different SU{N) groups. In fact, the 
fitted coefficient 03^3 agrees well with the estimate given in Eq. ( |3.2|) . This estimate 
does seem to favour the fit I. 

While the coefficient 04,3 only changes by a factor 2 in the case of SU{2) for 
both fits, the variation is about twice as large in the case of SU{3). However, it will 
become clear from the discussion in the following section that already our present 
estimates for 04^3 are rather restrictive for the effective mass of magnetic excitations 
in the plasma phase, since this latter is only sensitive to its cubic root. 



4 The O(g^) coefficient and its interpretation 



Combining ( p.7|) , (|3rT|) ,( |3.3| ) and ( |3.4| ) for case II we relate the lattice and continuum 



coefficients in the following way: 



a'. 



(2iV)3"4' 
1 



«4 = --^fc4,3 + cl(log(2iV) + i)). (4.1) 



Case I corresponds to setting C4 = 0. The values of 04,3 and C4 provide the input 
from the magnetic sector for the determination of the 0{g^, log g) terms in the 
weak coupling series of the free energy density. 
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The numerical values of the previous section lead to: 




(0.010 ±0.001), SU(2) 
(0.07 ±0.02), SU(3) 



(4.2) 



and 




(0.002 ±0.016), SU(2) 
(0.17 ±0.06), SU(3) 



(4.3) 



For the SU(2) group the estimated value for 04 in case II is compatible with zero 
within errors. Clearly, this means that the accuracy of our simulation in the SU(2) 
case is yet insufficient to become sensitive to logarithmic corrections. We note, 
however, that in the SU(2) as well as in the SU(3) case the inclusion of a possible 
logarithmic term in the fit at most increases the value of 04 by a factor of about 2.5. 

Let us try to get some feeling for the magnitude of 04 and its interpretation in 
terms of physical excitations. The 3-d theory represents the magnetic fluctuations 
of the high-T (3+1) -dimensional gauge theory. It is natural to expect that the 
0{g^) contribution to the free energy can be viewed also as resulting from a weakly 
interacting (almost ideal) gas of some pseudo-particles of mass uim ~ g^T and 
degeneracy Nd- This is quite analogous to the electric sector, where the 0{g^) 
contribution to the free energy density is just the contribution of a massive free gas 
of pseudo-particles with mass ~ gT. The well-known 3-d, 1-loop vacuum energy 
has the cut-off independent part: 



In the last equality the proportionality of uim to g^ has been exploited. 

The picture of an ideal gas of pseudo-particles cannot account for the logarithmic 
piece. If present at all this term would arise from higher order interactions of the 
effective degrees of freedom. The natural combination appearing in the argument 
of the logarithm is the ratio of the ultraviolet scale (the lattice constant) and a 
dynamically generated infrared scale ttiq^. This scale is proportional to g^: ma — 



£3 = 



SNomlj dniM 
127r dgf 



Npmli 



(4.4) 
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cgS'I- This fixes the separation of the terms proportional to and to (7^ log (7, 
necessary for the quantitative investigation of the non-logarithmic piece: 



nonpert 



C4 + + ln(2ArcG)) 



(2iV)s 

The comparison of ( [4 .41 ) to ( |2.7| ) via ( |3.1| ) leads to 



(4.5) 



127ra 



nonpert ' V 1/3 
4 1 „2 



""=1 ^) 

At least two simple cases can be put forward for the pseudo-particle excitations 
resulting from the high-T magnetic modes. With the ansatz given by Eq. ( [4.6|) 
they correspond to different choices for the degeneracy factor N^. The extreme 
cases clearly are to set equal to the number of gluonic degrees of freedom or 
to unity. The first case may be considered as being the analogous interpretation to 
the electric sector, while the second assumes that there do exist only color singlet 
magnetic modes above T^. 

i) Adjoint SU(N) multiplet of screened gluons 

Our previous determination of the magnetic screening mass of the gluons has 
been performed in Landau gauge 0. For reasons of comparability we continue to 
use this gauge. Then the magnetic mass and the infrared regularization mass are 
naturally identified. The number of degrees of freedom in this gauge is then given 
by Nf) = 2{N'^ ~ !)• The solution cg of the equation ( [4.6|) leads to the following 
mass estimates: 

(/) ^ 1(0.397 ±0.008) SU(2) 
"^giuon I (0.55 ± 0.04) SU(3) ^ ' 

and 



m 



{II) 

gluon 



(0.78 ±0.29)^1, SU(3). (4.^ 



This later value is still compatible within error bars with case I. Within the errors of 
the simulation the SU(2) result is perfectly compatible with the gluon mass calcula- 
tions presented in Ref. |Q. We also note that the ratio of the SU(2) and SU(3) mass 
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values for case I is very close to 2/3, which is expected from the 1-loop gap equation 
approach of since the 1-loop diagrams contribute a quantity proportional to N 
to the self energy. 
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ii) Scalar, SU(N) singlet glueballs {Nd=1) 

One arrives at the following prediction for the glueball mass: 

(/) ^1(0.721 ±0.014)^1, SU(2) 
ri-isingiet j ( 1 .38 ± 0. 1 1) SU(3) ^ ' 

and 

miZliet = (2.08 ± 0.67)^?|, SU{3). (4.10) 

For the determination of in case II one requires as additional input the infrared 
regularization scale. For the SU(3) case we used the measured magnetic screening 
mass, determined for SU(2) 0, and scaled it up by a factor 3/2, that is we used 
Cg ~ 0.6 in Eq. (|3p. 

These values can be partly compared with existing numerical results on 3-d 
glueballs. For the SU(2) theory the lowest glueball mass calculated in numerical 
simulations rrigiuebaii = 6.34(6)5^1 is clearly much larger than our estimate. Additional 
(higher) glueball states would make this discrepancy even more dramatic. Also the 
estimate of zero temperature glueball masses in four dimensions and estimates of 
finite temperature glueball screening masses in the (3+l)-d SU(3) gauge theory lead 
to larger values |jl5|, though the value found in case II within errors is on the edge 



of being compatible with the MC estimate. 

The above comparisons seem to present a rather strong evidence against an 
interpretation of the thermodynamics of magnetic fluctuations in terms of singlet 
excitations with a mass similar to a typical zero temperature glueball mass. This 
analysis suggests that the relevant thermodynamic degrees of freedom at high tem- 
peratures in the magnetic sector of the non-Abelian plasma are screened elementary 
gluons. 

5 Conclusions 

We have performed a first analysis of the non-perturbative contribution to the 0{g^) 
coefficient of the free energy of a finite temperature SU(N) gauge theory. We stress 
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that the analysis we have presented here is not complete. At the present stage 
we tried to investigate whether a lattice calculation can provide the necessary non- 
perturbative information that is needed to determine the complete 0{g^) term in 
the free energy of an SU{N) gauge theory. Our calculation shows that this is 
indeed feasible. Combined with a rather straightforward calculation of the expansion 
coefficient 03^3 within lattice perturbation theory and an analytic determination of 
C4 a much better determination of 04 can be achieved already on the basis of the 
numerical results presented here. The quality of the numerical input may also 
be further improved by exploring improved actions which reduce the size of the 
strongly cut-off dependent parts of observables like the expectation value of the 
Euclidean action. It is very encouraging, however, that already at the present stage 
we can distinguish between qualitatively different ideas on the nature of the high 
temperature magnetic excitations. 
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